Increasing identical particle entanglement by fuzzy measurements 
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We investigate the effects of fuzzy measurements on spin entanglement for identical particles, 
both fermions and bosons. We first consider an ideal measurement apparatus and define operators 
that detect the symmetry of the spatial and spin part of the density matrix as a function of particle 
distance. Then, moving on to realistic devices that can only detect the position of the particle to 
within a certain spread, it was surprisingly found that the entanglement between particles increases 
with the broadening of detection. 
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Entanglement is one of the most studied subjects in 
contemporary physics. This quantum mechanical prop- 
erty is the core of a whole new area in physics, mathe- 
matics and computation, known as quantum information 
theory. It has also been perceived as a useful resource for 
different technological applications 0. However much 
has still to be done in order to completely understand 
these non-classical correlations. One of the important 
aspects that still deserves further investigation concerns 
the actual observation of entanglement due to the indis- 
tinguishable character of the quantum world. 

Usually, when talking about entanglement, one tends 
to ignore the role of the measurement apparatus, always 
considering ideal situations. However, there is no such a 
thing as an ideal detector, and the detector bandwidth 
affects the measurement of entanglement. Furthermore, 
in the particular case of identical particles, it is still not 
clear how the symmetry of detection in external degrees 
of freedom affects the entanglement of the internal ones. 

In this paper we analyze these questions in two dif- 
ferent systems: a gas of non-interacting fermions at zero 
temperature and a photonic interferometer. In particu- 
lar, we analyze how the measurement of external degrees 
of freedom can affect the entanglement in internal de- 
grees of freedom in a system of identical particles. First, 
we show how the symmetry of detection affects the en- 
tanglement of two fermions detected in a non-interacting 
fermion gas at zero temperature. Then, we study how im- 
perfect detections affect the observation of entanglement 
both for the fermionic gas and also for an interferometer 
of photons. 

Recently Vcdral showed that the entanglement be- 
tween two particles in a non-interacting Fermi gas at 
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zero temperature decreases with increasing distance from 
each other 0] ■ Vedral also shows that there is a limit be- 
low which any two fermions extracted from the gas are 
certainly entangled. In particular, if both fermions are 
extracted at the same position, then the Pauli exclusion 
principle forces their spin to be maximally entangled in 
a typical antisymmetric Bell state. 

In order to calculate the upper limit for guaranteed 
entanglement, the author used the Peres-Horodecki cri- 
terion which in his paper, corresponds to f 2 > |, 
where f(r — r') = 3ji(kp\r — r'\)/(kp\r — r'\), \r — r'\ 
is the distance between the fermions, Pf = Kkp is the 
Fermi momentum and ji is a spherical Bessel function. 
In order to obtain this result, the author defines (apart 
from normalization) the spin density matrix of these two 
particles from the second order correlation function (see 
the Appendix for a justification) 

P l 3,w = ($o\H(r')*l(r)*Ar')* s (r)\$ ), (1) 

which corresponds to measuring these two fermions at 
positions r and r' (the subscriptions s,s', t, and t' refer 
to the spin of the particles). The initial state \4>q) is the 
ground state of the Fermi system 



n s , P bi( P ) io) , 



(2) 



and the detection operator at position r and spin s is 
given by, 

w d 3 P 
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pr bl(p). 
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A typical interpretation of this result is related to the 
spatial overlap of the wavefunctions of each fermion. 

Here, we present a clearer and more complete explana- 
tion to this behavior through the analysis of the symme- 
try of the position detection. To do so let us define new 
detection operators: 

UtAr,r') = (4a) 

n-,(ry) . (*.WAr')-*Ar)*.{r>)) ^ 
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The operator n+, (II SS ,), detects the antisymmetric 
(symmetric) spatial part and the symmetric (antisym- 



metric) spin part of fermion wavefunction. 

In terms of these new operators, Eq.Q) becomes: 



r 



|5>n 



(5) 



Note that written in this way, the two-fermion density 
matrix is the sum of four different terms, two of which 
contain only symmetric and antisymmetric spin detec- 
tors. The other two are the crossing terms, which vanish 
due to the exclusion principle (spin and position of two 
fermions cannot be both symmetric or antisymmetric). 
The remaining terms are: 

Psym = ^(*o|n+ t (r,r')II+ (r,r')|*o), (6a) 

Pasym = \ (*o| ^(r, r')^ , (r, /) |* ) • (6b) 

Here p sym takes into consideration only the symmetric 
spin function (therefore, it is related to the detection 
of the antisymmetric part of the spatial wavefunction) 
while Pasym contains only the antisymmetric spin func- 
tion (therefore related to the detection of the symmetric 
part of the spatial wavefunction). The density matrix 
can be rewritten as: 



P 



(2) 



Pasyjn H - Psym 

dpdp'^ { (1 + e i ( p - p 'K r - r ">) 



1 -1 
-1 1 



+ (1 — e i (p-p')( r - r ')\ 



( 2 
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1 1 
1 1 



(I) 



First, note that, as expected, for r = r' . the antisym- 
metric spatial function goes to zero, and the spin wave- 
function has to be antisymmetric (first term in Eq.Q). 
For r — r' ^ both parts contribute. Note, however, 
that the symmetric spin density matrix can be viewed as 
an equal weight mixture of the three triplet components, 
and has no entanglement at all (by the Peres-Horodecki 
criterion). The spin state in Eq.Q represents a con- 
vex combination of singlet state and the equal mixture 
of triplet states. It will have entanglement iff the frac- 
tion of singlet is sufficiently larger than the fraction of 
triplets, in a similar context to the calculation of relative 
robustness In the limit r — r' 3> l/kp, as the inte- 
grals on Eq.Q are performed over momenta below the 
(momentum equivalent of) Fermi surface the momentum 
dependent terms oscillate too fast, and average to zero. 
Therefore, the spin density matrix becomes just the iden- 
tity (StsSt's')- This behavior can be seen as a smooth 



transition from a quantum statistics (Fermi-Dirac) to a 
classical one (Maxwell-Boltzmann). 

Another way of analyzing Eq.Q is given in Ref.Q, 
where the authors describe a"' MS (X Werner state. Here, 
it becomes clear how the symmetry of the spatial detec- 
tion is influencing the behavior of the spins, and also why 
there is no entanglement from \r — r'\ ~ l/kp on. 

An interesting question arises when considering non- 
punctual (non-ideal) position detections, i.e. more real- 
istic apparatus that detect the position of those fermions 
with some uncertain. Instead of Eq.@, the detection 
operator should be written as the general field operator: 



* a (r) = jj e lpr "D(r - r")b s (p)dr" dp. (8) 

The perfect position measurement situation is the partic- 
ular case of Eq. JSJ) corresponding to D(r—r") = 5(r—r"). 
However, if D(r — r") has some position uncertainty, like, 
for example, if it is described by a gaussian with spread 
o-, 



D(r - r") 



1 



2na 



0) 



then, the field operators become 



tf,(r) = -^=- [[ e ipr "e 'aX'" b s (p)dr"dp. (10) 



27TCT 

These operators, when substituted back in Eq.QJ, give: 
Ps S ',w = 5 st 5 s > t >f(d,a)+5 st >5 a > t g(d,a), (11a) 
with d = \r — r'\, 

'erfjapfy 2 



g(d, a) 



f(d,a) 



1 -d 



(lib) 



id , 



id 



3^ {erf (op/ - — ) - erf(-— )} 



2a' 

x {crf(cr/ + — ) -erf(— )}, 



2a' 



2a' "2a' 
and erf(x) is the "error function", defined via: 



erf(x) 



exp(—t 2 )dt. 



(11c) 
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In order to make clearer the behavior of entanglement 
in relation to changes in d and a we can compute the 
negativity (M) @ of the state p^(d,a,pf). This en- 
tanglement quantifier was chosen because it is an en- 
tanglement monotone (i.e.: satisfies all the required 
features of a good entanglement quantifier), and also it 
is easy to calculate. Af(p) is equal to twice the absolute 
value of the negative eigenvalue of the partial transpose 
of p. p TA . N(p) measures by how much p TA fails to be 
positive semi- definite 0| In our case M(p) reads as 



Af( P ) = 2 M ( 



/-2g > 
4/-2<?' 



(12) 



where 



-x , x < 
, x > 
Fig. 2] for some values of a. 



This function is plotted in 
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FIG. 1: (Color online) M{p) vs. d for p; 
(grey), a = 2 (magenta) and a — 4 (black). 
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Note that, for imperfect position detection, the entan- 
glement decreases as the detectors become apart from 
each other, but increases if the spread in the detection 
becomes larger. The fact that inaccuracy in the detec- 
tion increases entanglement seems surprising. However 
it has to be noted that as our knowledge in position gets 
worst, our knowledge in momentum gets better. In the 
limit of infinite spread, both detectors become perfect 
momentum detectors (centered at p = 0, see Eq.Q), 
which means again that their spin wavefunction should 
be totally antisymmetrized, hence they are found in the 
antisymmetric Bell state. It is important to stress that 
Eq. (JSJ describes a coherent combination of localized field 
operators instead of a statistical average of them. That is 
the reason for the infinite spread limit be a momentum- 



localized detector instead of just a vague "there is a par- 
ticle somewhere" . 

A similar situation can be found in a bosonic sys- 
tem calculating the entanglement in polarization between 
photons placed in an interferometer using polarization 
beam splitters (PBS). In this analysis the polarization 
and frequency of the photons play the role of the spin 
and momentum of the fermionic case. This kind of in- 
tcrferometry has been extensively studied in practical ex- 
periments Q and thus, in principle, the results presented 
here could be tested experimentally. 

First we write the modes after the PBS if the modes 
before it are ai j£ and a 2 ,e (see Fig|2j): 

£ 3 ,e(w) = T e a he {w) + iR t a 2 ^[uj)e^ T , (13a) 
E i<t (u) = iR eai . e (uj)+T*a 2 . e (uj)e lu,T , (13b) 

with |T e | 2 and |i? e | 2 being the transmissivity and reflec- 
tivity of the e polarization in the PBS, r = t\ — t 2 , where 
Ti and T2 are the propagation times from the crystal to 
the detectors via path 1 and path 2 respectively. 
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FIG. 2: (Color online) Hong-Ou-Mandel type of interferome- 
ter. 

The two particle density matrix is again defined 
through the second-order correlation function: 

Pee>,ss> = JJ dui 1 dui 2 D(uJ 1 ,U}2) (14) 
X (0 O | £l,>l)£| )S ,(cJ2)£4,e^2)£3,e(wi) |<fo) . 

The state 

\<M - ^ 10) (15) 
is prepared by the two-photon creation operator 

^ = ^ a\ k (i' 1 )al k ,(iy2)gkM'(^i,L'2)d^idv2, (16) 
k,k' 

with frequency distribution gk,k\v\ 1 V2) which is, in the 
more general case, polarization dependent. Note that 
this calculation is more general than a typical Hong-Ou- 
Mandel interferometer, since one is free to choose initial 
states by gk,k'( u ii v ?)- I n a w &y similar to Eq.JSJ), detec- 
tion is represented by a function D(uj\,lo2). 
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Using the commutation relation for bosons and eval- k, k' we find: 
uating the integrals in dv\ , dv-i and the summations in 



Pu',s' s = JJ ' {T s >T*T*,T € g* s ,(u>l,U2)gee>{ui,V2) + R s >R* s R*,Rtg* s , s (u)2,Vl)ge'e(u2,Ul) (17) 

The PBS is an optical element which allows transmission in only one polarization (say vertical), and reflection in 
the other one (horizontal). The coefficients of transmissivity and reflectivity can be written as Rh = T„ = 1 and 
Th = R v = 0. With this, and assuming a symmetric state in polarization (g ss i = g s ' s ), the unique non-vanishing 
coefficients arc: 



Phh.vv 



\g t e{<-oi, w 2 )| 2 -D(^i, u 2 )du>iduj2, e = h,v, 
g* v (ui,ui2)ghh(u2,uJi)e l ( u ' 1 ~ UJ2) D(uJi,U2)dujiduj2 = p% Vt hh- 
I 



(18a) 
(18b) 



In order to carry on the calculations, let us assume 
particular distributions for the functions g(uji,oj2) and 
D{uji,u>2)- Since we are more interested in the detection 
properties, let us first assume that g(u>\,uj2) is constant. 
Then, we assume independent gaussian detectors, which 
means D(oji,uj2) = D(uj\)D(lo2) and 



D(w) 



1 



2na 



(19) 



Now, we can calculate the polarization density matrix, 
obtaining, up to normalization, 



Phh,hh 

and 

Phh,vv 



,j,w = JJ D(ojx)D(uj2)du)i<L)2 



Pvv,hh 



1 r < w 
w J exp(_ ^ 

I exp(-^-^- - iuJ2T)dbJ2- 



(20) 
(21) 



Evaluating these Fourier transforms of Gaussian func- 
tions we obtain 



Phh,vv Pvv.hh 



(22) 



The (normalized) polarization density matrix then is 
given by 



P= 2 



I I 0-/ 



V -/ o o i 



(23) 



where / = exp(— <t 2 t 2 ). The partial transposition of p 
has always one negative eigenvalue equal to — ^. So the 



Negativity of p is Af(p) = f = exp(— <t 2 t 2 ). This result 
shows that as the detector becomes broader in frequencies 
or the time delay between the two arms of the interferom- 
eter becomes larger, the entanglement between detected 
photon pairs decreases. One should note that time delay 
gives rise to which path information, while broadening in 
frequencies opposes to such labeling. 

Again we can see the state p as a combination of Bell 
states. Interestingly, in the intcrferomctric process de- 
scribed here the separation is between the states |$+) 
and |$— ), in the following way: 



P = 



(1-/) 



(1 + /) 



(24) 



and no more between symmetric and antisymmetric 
states as in Equation J7J. 

In conclusion it was shown that the measurement ap- 
paratus plays a central role in the entanglement of iden- 
tical particles. That entanglement increases because of 
broadening in detection can sound weird at first. How- 
ever it is important to say that the detections discussed 
here are done in a coherent way (as in most of real cases). 
If we have modeled the imperfect detectors as a cluster 
of perfect detectors things would be different. Further- 
more the relation between identical particle entanglement 
(i.e.: the correlations due exclusively to indistinguisha- 
bility) and the entanglement generated by direct interac- 
tion among the particles is not completely understood. 
Despite of that the framework presented here can be, a 
priori, also applied for interacting particle systems and 
then could make this relation less obscure. 
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APPENDIX A: SECOND ORDER CORRELATION 
FUNCTIONS AS QUANTUM STATES 

In this appendix we discuss why we consider the sec- 
ond order normalized correlation functions and 
as two-qubit quantum states. Mathematically speaking, 
any trace 1 positive operator is a state on the space state 
in which it is defined. Physically, however, this seems a 
poor justification. Moreover, why one can consider this 
four-dimensional linear space as a two-qubit system and 
obtain conclusions about its entanglement? 

To answer this questions we need to consider an in- 
terplay of two aspects of quantum theory: photodetec- 
tion and post-selection. In a first quantized notation, 
if the electromagnetic field state is described by p and 
we make a one-photon detection described by the anni- 
hilation operator a, the field state after such measure- 
ment will be given by p' = apaJ , properly renormal- 
ized (i.e.: the conditional state). As photodetection is 
a destructive measurement, after the click, p' describes 
the state of the remaining photons. If, however, we 
have a non- destructive one-photon detector, this post- 
selected "detected photon" would be described by the 



pure Fock state |1) of the mode annihilated by a. But 
what if we just have a non-destructive one-photon de- 
tector similar to the one above, except by the fact that 
it clicks for any one-photon superposition of two specific 
modes (e.g.: two polarizations of one and the same spa- 
tial mode). This situation is analogous to a degenerate 
spectrum projective measurement. The state of the post- 
sclcctcd photon will preserve the coherence which it pos- 
sibly had between the two modes, being described by a 
linear combination like a |1, 0) + f3 |0, 1), or, more gener- 
ally, by a density operator p = PpP, with P the projector 
P = |1, 0) (1, 0| + |0, 1) (1, 0|, and again properly normal- 
ized. Clearly, the state of this two-level system describes 
the unresolved mode combination (in the example, the 
polarization degree of freedom). 

If we detect two spatially distinguishable photons, but 
each of them with unresolved polarizations, the first 
quantized notation becomes cumbersome, but the ideas 
involved are the same, and the remaining two-qubit state 
describes the polarization degree of freedom of this two 
distinguished photons. 

The same reasoning, with two fermions spatially lo- 
calized ar r and r' , but with unresolved spin, in sec- 
ond quantized notation, give rise to Eq.(Q), that then 
describes the spin state of the two localized fermions [Tol|. 
In this sense, (normalized) second order correlation func- 
tions can be considered as quantum states, and its en- 
tanglement properties can be studied, as is done in this 
paper. 
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